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We study manipulation of entanglement between two identical networks of quantum mechanical 
particles. Firstly, we reduce the problem of entanglement transfer to the problem of quantum state 
transfer. Then, we consider entanglement concentration and purification based on free dynamics on 
the networks and local measurements on the vertices. By introducing an appropriate measure of 
efficiency, we characterize the performance of the protocol. We give evidence that such a measure 
does not depend on the network topology, and we estimate the contribution given by the number of 
entangled pairs initially shared by the two networks. 
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I. INTRODUCTION 

Networks of quantum mechanical particles play a fun- 
damental role in extending quantum information process- 
ing to multipartite settings and in foreseeing the realiza- 
tion of distributed protocols for multi- users devices [T] . 

In the present work we study manipulation of entan- 
glement between quantum networks and attempt to iso- 
late useful graph-theoretic properties in relation to the 
dynamics of entanglement. Specifically, we analyze the 
problem of entanglement transfer through networks and 
then we relate its setting to the fidelity problem associ- 
ated with state transfer. 

Our goal is to study entanglement concentration and 
purification, when the evolution of the system involves 
two networks of spin-half particles with Xy-type inter- 
action. Entanglement concentration (and, more gener- 
ally, purification) is among the most relevant tasks in 
quantum information processing [2J. It is a process ex- 
tracting strongly entangled pairs out of initially weakly 
entangled ones, only by mean of local operations and 
classical communication. Most of the known protocols 
involve control- not operations. In a similar way, mul- 
tipartite measurements and classical communication are 
used to propagate entanglement in certain type of quan- 
tum networks [3, Entanglement concentration has 
also been obtained by exploiting the effects of quantum 
statistics of indistinguishable particles [1]. A streaming 
(sequential) protocol for universalj;ntanglement concen- 
tration has recently appeared in [6], where a number of 
input systems are processed in sequence, and perfect en- 
tangled pairs are obtained at the end of the protocol. 

Following 0], we simply exploit free dynamics on the 
networks and local measurements on the sites. Somehow 
contrarily to a claim proposed in Q, we show that this 



method works reasonably well also in the case of XY- 
type interaction. We consider generic network topologies 
other than linear spin chains and we introduce an effi- 
ciency measure to characterize the performance of the 
protocol. It will turn out that, at least in the cases con- 
sidered, the structure of the networks only affects the 
time at which the maximum efficiency is attained. 

The topic of this paper fits into a growing literature. 
State and energy transfer in spin systems are now areas 
of extensive study ^ (see also 9]). Protocols for cre- 
ating entanglement between distant sites in a quantum 
network have been designed in 10], where the network 
considered is a lattices, whose vertices share pure, non- 
maximally entangled pairs of qubits. The problem of 
preparing maximally entangled states between different 
sites in a distributed quantum system is related to classi- 
cal percolation in statistical mechanics [ll|. In [l^], the 
study of entanglement percolation is extended to struc- 
tures with precise statistical properties but unknown lo- 
cal topology. In [Tsj . entanglement percolation is used 
for long-distance singlet generation in networks with ad- 
jacent vertices connected by partially entangled bipartite 
mixed states. 

The structure of the paper is as follows. In Section II, 
we introduce the considered networks model relying on 
the basic XF-type interaction. In Section III, we give 
rigorous results for the problem of entanglement transfer 
in quantum networks by considering the single excitation 
model. The quality of the entanglement transmission is 
then related to the fidelity achieved between the con- 
sidered nodes, therefore extending the cases studied in 
Q. In Section IV, we discuss how to enhance the en- 
tanglement of a pure states through free dynamics and 
local measurements on the nodes. We introduce an effi- 
ciency measure for our protocols, by considering all possi- 
ble measurement outcomes and the corresponding prob- 
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abilities. Finally, in Section V, we show how the protocol 
works in the case of mixed entangled states by consider- 
ing Werner mixtures Likewise Section IV, the ef- 
ficiency of the procedure will be related to the minimal 
number of initially entangled pairs necessary to success- 
fully complete the protocol. Conclusions are drawn in 
Section VI, where we summarize the results and mention 
potential physical realizations. 



II. NETWORKS MODEL 

Let G = (y, E) be a simple undirected graph (that is, 
without loops or parallel edges), with set of vertices V{G) 
and set of edges E{G). The adjacency matrix of G is de- 
noted by A{G) and defined by [A{G)]^J = 1, if G E{G)- 
[A(G)]i,j = if ij ^ E{G). The adjacency matrix is a 
useful tool to describe a network of n spin-half quantum 
particles. The particles are usually associated to the ver- 
tices of G, while the edges of G represent their allowed 
couplings. If one considers the XY interaction model, 
then G E{G) means that the particles i and j in- 

teract by the Hamiltonian [H{G)]ij — — {XiXj + YiYj), 
where Xi, Yi are the Pauli operators of the i-th parti- 
cle (here we consider unit negative coupling constant). 
Thus, the Hamiltonian of the whole network reads 



(1) 



1 ®n 



and it acts on the Hilbert space H. = 

In the following we will consider two networks with 
underlying graph Ga and Gb having, for the sake of 
simplicity, an equal number of nodes n. Then, if we apply 
XY model we have two Hamiltonians 

1 " 

^(Ga) = --Y,[A{GA)h{XA^XA, + YaYa,). 
1 " 

H(Gb) = "2 Y.iAiGB)hiXB,XB, + YbYb,). (2) 



III. ENTANGLEMENT TRANSFER 

Let A and B be two networks sharing a pair of qubits 
in a maximally entangled state 



^(|0)|1) + |1)|0)). 



(3) 



Suppose that the entangled qubits are in the i-th site of 
the network A and in the fc-th site of the network B. We 
assume that all the remaining qubits are in the ground 
state. Then the state for the joint system A + B will be 



1 



V2 



mA\k)B + \l)A\0)B), 



(4) 



where \i)A denotes the state |0)i|0)2 ■ ■ ■ ■ ■ . |0)„_i|0)„ 
for the network A; a similar notation is used for the net- 
work B. 

The evolution of the whole system will be governed by 
the Hamiltonians ([2]). Hence, the global evolution of the 
two networks will lead to the state 



H(GB)t 



V2 



\k)i 



i)a\0)i 



V2 



(|0)A|</'fe)B + |V^^)A|0)s), 



(5) 



where = e'"^'^^'^^\i) a and \Lpk)B = e'"^-^B^^\k) b- 

To see how entanglement propagates from the sites 
i and k, we evaluate the concurrence [isj between two 
generic sites j and I (different from i, k). To this end, we 
have to consider the following density matrix pab'- 

PAB = iTry{|0)A(0|}$5Tr^|^fc)B(¥'fc|} 
+ iTry{|0)A(V'JI®Tr^l¥'fc)B(0|} 



-Try{|^,)^(0|}®Tr^|0)B(^fc|} 
iTry|^,)^(^,|(»Tr^|0)s(0|}, 



(6) 



where Try denotes the trace overall qubits but the j-th 
one. 

For a single pure excitation of a spin-half network, the 
wave functions can be expanded as 



= ai|l) 



+ an\n), 



or 



1^) =/3i|l) + ... + /3„|n). 



(7) 



(8) 



In this way it is possible to express the density matrix 
Pab as a function of the complex coefficients a and /3: 



PAB 







V 



\ 

lAP a*/3, 
a,/3; |a,f 


/ 



For the single excitation model including both net- 
works with Hamiltonians ([2]), the pab Pab matrix is 



PabPab 




2|a,flAI^ 



|/3,|2(a*/3, + «,/?;) 



2|a, 



V 















Here, accordingly to [l5| 

PAB = {Ya ® YB)p*ABiYA ®Yb), 



I 

(9) 
(10) 
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where p\g is the complex conjugate of pab [ISj . 

The concurrence between two qubits j and I is defined 

as fii 



C-j^i{p) = max ( 0, 



(11) 



where Ai, A2, A3, and A4 are the eigenvalues of pab Pab 
in the decreasing order. 
Explicitly, 



0,1 



(12) 



It turns out that the concurrence is equal to the product 
of the fidelity reached during the transmission of a single 
excitation in each network. As it is expected, if there is 
perfect state transfer between particle i and particle j on 
A as well as between particles k and I on B, entanglement 
is perfectly transferred from the pair i, fc to the pair j, I. 



IV. ENTANGLEMENT CONCENTRATION 

We now consider two identical networks sharing pairs 
of qubits in a pure, non- maximally entangled state and 
look for a possibility of enhancing the amount of such 
entanglement. 



A. Two qubits networks 

In the simplest configuration we may consider two sites 
on each network connected by an edge. These networks 
are simply two graphs consisting of a single edge. Then, 
the resulting Hamiltonians Ha and Hb correspond to the 
adjacency matrices 



A{Ga,b;P2) = 



1 

1 



(13) 



where P„ denotes an open chain of length n — 1. Let the 
initial state of qubits Ai and Bi be 



IV') = cos((?)|0)aJO)b, +sin(0)|l)^Jl) 



(14) 



and the initial state of qubits A2 and B2 be |0)yi2 |0)_b2 ■ 
The two networks will evolve independently, so that the 
whole state 

1*) =cos(0)|O)aJO)bJO)a,|O)b, (15) 
+ sin(^)|l)^Jl)BjO)^jO)B, 

will have the following general form at time t: 

|M'(i)) =cos(0)|O)^JO)^,|O)bJO)b, 

+ sin(0)cos2(t)|l)^jO)^Jl)BjO)B, 

+ sm{6) cos(i) sin(t)|l)^, \0)a, |0)bJ1)b2 

+ sin(0) cos(t) sin(t) |0) 1 1) a J 1)b, |0)b, 

+ sin{0) sin2 (i) 1 0) A J 1) A. |0) B J 1) . (16) 



Actually, given the evolution, also the qubits A2 and 
B2 become entangled with Ai and Bi. Local measure- 
ments on A2 and B2 may concentrate the entanglement 
available at time t on qubits Ai and Bi, possibly in- 
creasing the initial amount. Suppose that we locally 
measure Za2 ^ Zb2- For outcomes 00, the state in Eq. 
([T6| will be projected onto |0)a2|0)_b2- Hence, the state 
a = |*(i))(*(t)| changes into 



cr ^ PA^B^ «) |0)a2(0| «) |0)j 



where 



and 



It results 



PAiBi 



PAiBi 



TrAiBi{pAiBi}' 



PA^B^ ^A, {0\b.AOW\0)bM 



A2- 



(17) 



(18) 



(19) 



PA,B, =A2 (0|B2(0||«'(i))(*(t)||0)B2|0)A2 (20) 

= (cos(0)|O)aJO)b, -Ksin(0)cos2(i)|l)^Jl)sJ 
{cos{9)a, (0|si (0| + sm{0) cos^{t)A, {1\b, (1|). 

The probability of this transformation is 

Pooit) = Tta,bApa,bA = cos2((?) +sin2(0) cos^t). 

(21) 

Notice that the state dH]) with (HH) and is quite 
generally entangled with a concurrence possibly higher 
then the initial state (ITC)) . On the contrary, it is easy 
to check that the state resulting from the projection of 
Eq.dUl) onto |0)a2|1)b2) or |1)aJ0)b2, or |1)a2|1)s2 (cor- 
responding to outcomes 01, 10, 11 respectively) is sepa- 
rable, hence with zero concurrence. 

Since the enhancement of the initial entanglement is 
determined by the measurement outcomes, we have to 
account for their respective probabilities in order to eval- 
uate the performance of the protocol. Then, we define the 
efficiency of the protocol as the average of the increments 
of entanglement weighted by the respective probabilities, 
that is, 



E:=J2Po^oC. 



(22) 



Here the index o runs over all possible measurement out- 
comes, each occurring with probability Pq. Moreover 



(Co - C), if Co > C; 
0, if Co < C. 



(23) 



Here C is the concurrence of pa^Bx and Co the concur- 
rence of p'a^Bi conditioned to the outcome o. 

In our case the concurrence for the outcome 00 reads 



Can — 



2sin(6l) cos(6') cos^{t) 
cos2(6l) -|-sin2(e') cos4(t)' 



(24) 
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thus the efficiency simply results 

E{e,t) ^ PooAooC (25) 
= 2 sin(e') cos(6l) (cos^ (t) - cos^ (9) - sin^ (6) cos^ (t) ) . 

Now taking ^"^g^'*-* = we get the optimal time (maxi- 



mizing ((25|) ) as 



cos^{t) = 



1 



2sin^(6') 



(26) 



By substituting it back to Eq.(|25|) we arrive at the max- 
imum efficiency written as solely function of the variable 



2 sm(6') 



(27) 



The behavior of this quantity is depicted in Fig. (IT|). 

Let us now consider the case of two initially entangled 
pairs. Let Eq. ([T4l) be the initial state of qubits Ai and 
Bi as well as qubits A2 and i?2- Write VA^Bi = \'ip){'4'\: 
then the initial state of the networks simply results 



VAiBi (Xi VA2B2, 



(28) 



and it evolves into 



_ „~t.HAt„~cHBt 



a = e — (va.b, va^b^) e'"^'e'"^\ (29) 

Now let us consider local measurements of the Pauli 
observables Za2'SiZb2, i-C-, only one out of the two initial 
pairs is locally measured. 

For outcome 0262 (with 02, &2 G {0, 1}) the above state 
will be projected onto 1 02)^2 1 ^2)^2 • Hence, 



cr i-^ PAiBi 102)^2(02! ® \b2)B2{b2 



where 



and 



PAiBi 



PAiBi 



TrAiSi{pAii3i}' 



PAiBi = {a2b2\o-\b2a2) 



(30) 



(31) 



(32) 



Outcome 02^2 = 00 occurs at time t with probability 



Poo{t) = TiaibApaibA 



(33) 



= - cos^(6')(2 cos^{e) - sm'^{9) - sin^ie) cos(4i)), 

and the corresponding conditional state results (in matrix 
form) 



PA^B^t) = 
/ cos^(6l) 



(34) 

cos3(6')sin(6')cos(2t) \ 
' 



cos3(6')sin(6')cos(2t) cos2(6') sin2(6') cos2(2t) 



Its concurrence reads 

4sin(6') cos(6') cos(2t) 



Cooit) 



(35) 



2cos2(6l) - sm^{9) - sin^{9) cos(4t) ' 
Outcomes 02^2 = 11 occurs at time t with probability 



Pll{t) = TrA^BiiPAiBi} 



(36) 



1 



= ^ sin2(6l)(cos2(6l) -I- -2 sin^{9) + cos^{0) cos(4i)), 

and the corresponding conditional state results (in matrix 
form) 



PAibA^) 



(37) 



/ cos^ (61) sin^ (6*) cos2(2t) cos(6')3 sin(6') cos(2t) \ 

' ' 



cos(6l)sin^(6')cos(2t) sin'' (6*) 



v 

Its concurrence reads 

4sin(6l)cos(6') cos(2t) 



Cii{t) 



cos^ (61) - 2 sin^ (61) + cos2 (9) cos(4i) 



(38) 



Since the other outcomes (0262 = 01, 10) give separable 
conditional states, the efficiency in Eq.(|22|) becomes 



E{9,t) = cos^(6i)sin'*(6') 

X (cos2(6l) - 2 sm^{9) - 2 cos(2t) + cos^{9) cos(4i)) 

X (2 cos^{9) - sin2(6l) - 2 cos(2i) - sin2(6i) cos(4t)) 

- 2cos^(6l)sin(6')cos(4t). (39) 



Now, by taking "^g^'^' = we get the optimal time (max- 



dE(e,t) _ 
dt — 

imizing (|39)) ) as i = tt. By substituting it back to Eq. 
P9p we arrive at the maximum efficiency written as solely 
function of the variable 9 



E{9) = 2 cos"* (61) sin(6l) (8 sin^ (61) - cos(6')) 



(40) 



The behavior of this quantity is depicted in Fig. ([T]) for 
7r/4 < 6* < 7r/2. 



B. Networks with more than two qubits 

In the previous subsection we have seen that having ini- 
tially a single entangled pure pair suffices to achieve en- 
tanglement concentration. However, looking at Fig. (fl]). 
we can conclude that having initially two entangled pure 
pairs give to the protocol a better performance. From 
now on we are going to consider a number of initially 
entangled pure pairs equal to the number of vertices in 
the network. For n vertices, consider the initial state of 
qubits Ai and Bi to be (fT4|) . Writing VA^Bi = 
the initial state of the system is then 



VAiBi (X VA2B2 



(41) 
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The time-evolution gives 



(42) 



Now we measure with respect to the Pauh observables 
{Za2 «) . . . (8) ZaJ <E) {Zb2 <S)-.-(E) ZbJ, i.e., n - 1 out of 
initially n pairs are locally subjected to measurement. 

When we obtain 02 . . . a„b2 where each variable 

is either zero or one, the above state will be projected 
onto [02)^2 • • • |a«)A„ 1^2)^2 • ■ • l^n)B„- Then 



o- PAiBi fXi 102)^2(02! (8) ... (8) \a„)A^{an\ 
(8 162)52(^21 (8... (8 |6„)b„(?'«|, 



where 



PAiBi 



TrAii3i{pAii3i}' 



and 



PAiBi 



(02 . 



■ anb2 . . . bn\cr\b 71 . . . 020-71 . . . (3-2 ) 



(43) 



(44) 



(45) 



By taking into account all 4"~^ possible outcomes we can 
then calculate the efhciency given in Eq. (|22p . 

Let us now consider n = 3 and the two networks in an 
open chain configuration. Then, the Hamiltonians Ha 
and Hb come from the following adjacency matrices 



1 

^(Ga,b;P3) - I 1 1 
1 



(46) 



The efficiency is maximized by t = kn/\/2, fc e N giving 

E{e) =:5cos2(6')sin2(0)(l-2cos2(6')), (47) 

where n/A < 9 < tt/2. 

The same efhciency can be obtained by considering the 
dynamics resulting from the fully connected three qubits 
networks. The related Hamiltonians Ha and Hb are 



A{Ga,b;K3) = 




(48) 



where is the fully connected graph with n vertices. 
However, in this case, the expression (|47p is obtained at 
time t = kn, k G N. The quantity (|47|) is represented in 
Fig. (P. 

Practically, the change of the underlying graph associ- 
ated to the network only affects the time of the dynamics, 
but it does not affect the expression of the maximum effi- 
ciency. We conjecture (supported by numerical evidence) 
that this is also true for networks having more than three 
nodes, provided that there are no isolated nodes. In fact, 
the protocol involves entangled pairs shared by single 
nodes of each network. Since we consider a number of 
initially entangled pairs equal to the number of nodes in 
a network, two isolated nodes that share an entangled 



pair would constitute a separate system. Then, the max- 
imum efficiency only depends on the number of nodes 
and so on the total number of initially entangled pairs. 

Let us now consider n = 4 and suppose that the two 
networks A and B are both in a closed chain configu- 
ration (a cycle). Then, the Hamiltonians Ha and Hb 
result from the following adjacency matrices 




(49) 



where C„ is the closed chain of length n—l. In this case 
the maximum efficiency is obtained for t = kir, k G N 
and it reads 

E{e) = Y^cos(6') sin3(6i)(8cos(46')-3cos(66')-29cos(26') + 

(50) 

where 7r/4 < 6 < 7r/2. Such efficiency is attained at 
t = tt/2 + kn, k e N, in the case of fully connected 
configurations. 

The efficiency ((50)) is represented in Fig. ([1]). 

Let us finally consider the case of n = 5. We chose 
Hamiltonians Ha and Hb coming from the following ad- 
jacency matrices 



A{Ga.b;C5) 



/ 1 I \ 

10100' 

10 10 
10 1 

V 1 1 0/ 



(51) 



The maximum efficiency results 



1 



E{e) = -cos(6l)sin3(6')(cos(46l)-cos(66l)-13cos(20) + l), 

^ (52) 
where 7r/4 < 9 < tt/2. This is obtained for t = kTT, k eN. 
The efficiency ([52l) is represented in Fig[TJ 



V. ENTANGLEMENT PURIFICATION 

We shall consider pairs of qubits in networks A and B 
initially in a mixed entangled state. We then determine 
conditions for enhancing such entanglement. Specifically, 
we consider the initial state of qubits AiBi to be the 
Werner mixture '14'I written as 



WA.Bi 



A,B, 



1-/ 



/AiB: 



where 



^+)A,B,(*+| + |*")A,i3.(*"|) , (53) 

|0)aJO)b. ±|1)aJ1)b. 



V2 

\0)aM)b,±\1)aMb, 
V2 



(54) 
(55) 
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FIG. 1: Maximum efficiency E{9) as a function of 9 for differ- 
ent networks. Curves from bottom up refer to networks with 
number of nodes and number of entangled pairs equal to (2, 1) , 
(2,2), (3,3), (4,4), (5,5) respectively. 



are the Bell states and / is the fidelity between w and 
|$+) (notice that 0.25 </<!). The initial state of the 
networks is 

Wi2...n = WAiBi ® WA2B2 ® • ■ • (8) Wa^B^ ■ (56) 

By using the same procedure shown for entanglement 
concentration, we can calculate the efficiency of the pro- 
tocol. Networks with one or two qubits lead to a zero 
efficiency. In fact, the concurrence obtained after the 
evolution of the initial state in Eq. ([TE| is less than the 
initial one and it is identical for every outcomes. Actually 
for all outcomes o we have 

Co - C = (1 + 4/-32/2)/36 0. (57) 

The efficiency for three qubits becomes greater than 
zero and its maximum is represented in Fig. Such 
maximum is obtained when t — kn/\/2, fc G N, in the 
open chain configuration (from (j46p and when t — kn, 
fc € N, in the fully connected configuration (from (|^5)) ). 

Likewise the case of entanglement concentration, only 
the time t at which maximum efficiency is attained is 
affected by the choice of the networks' topology (not the 
maximum efficiency itself), provided that there are no 
isolated nodes. 

The maximum efficiency for n — 4 and n — 5 can be 
obtained from Eqs. and (IFTI) . respectively. Their 

behavior is represented in Fig. ([2]). 



VI. CONCLUDING REMARKS 

In conclusion, we have considered entanglement ma- 
nipulation between two networks of spin-half particles 
with XY-type interaction. We have related entangle- 
ment transmission to the fidelity achieved in state trans- 
fer when adopting the single excitation model. 
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FIG. 2: Maximum efficiency E{f) as a function of / for dif- 
ferent networks. Curves from bottom up refer to networks 
with number of nodes and number of entangled pairs equal to 
(3,3), (4,4), (5,5) respectively. 



We have shown how to increase the entanglement of 
pure and mixed states via free dynamics and local mea- 
surements on the nodes. The performance of the protocol 
are characterized by means of an efficiency measure. The 
performance appears to be independent of the networks' 
topology provided that there are no isolated nodes. The 
structure of the network only affects the time at which 
the maximum efficiency is attained. Moreover, we have 
determined the minimal number of initial entangled pairs 
shared between the two networks and needed to success- 
fully complete the protocol {i.e., with nonzero efficiency). 
This number is 1 and 3 for entanglement concentration 
and purification, respectively. Furthermore, the results 
shown in Figs. Q and ([2]) suggest that a (quasi) con- 
stant increase for the maxima of the efficiency is obtained 
with a corresponding increase for the number of entan- 
gled pairs. This can be estimated as about 0.15, for the 
case of entanglement concentration, and about 0.04, for 
the case of purification for each additional pair. This 
trend suggests that a (quasi) deterministic protocol can 
be achieved, at least for specific values of 9 or /, with a 
limited number of initially entangled pairs. Our protocol 
is straightforward and only based on the free dynamics 
occurring in the two networks. It may be an interesting 
venue to consider more general settings, as, for example, 
synchronous dynamics in different networks. 

Notice that already in the single excitation model the 
XYZ (Heisenberg) interaction would give a different net- 
works' dynamics. Actually a phase factor would be intro- 
duced in front of the state's components containing one 
excitation per network. However such a phase factor is 
not influent to our protocol as long as we consider mea- 
surements on the computational basis (Z eigenvectors). 
Thus our conclusions can be applied to XY model as well 
as to Heisenberg model. 

Finally, it is worth remarking that the proposed 
schemes for entanglement manipulation might be experi- 
mentally verifled with networks containing few spins-half 
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systems. Candidate physical systems would be semicon- 
ductor quantum dots [l^|, superconducting qubits [17] or 
atoms trapped in optical lattices [l8{ . 
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